FOUR- AND FIVE-LAYER SILICON-CLAD 
DIELECTRIC WAVEGUIDES'^* 


Glen McWright and T. E. Batchman 
University of Virginia 
Charlottesville, Virginia 


ABSTRACT 


Computer modeling studies on four-layer silicon-clad planar dielectric wave- 
guides indicate that the attenuation (a) and mode index (3/K) behave as exponentially 
damped sinusoids as the silicon thickness is increased. The observed effect can be 
explained as a periodic coupling between the guided modes of the lossless structure 
and the lossy modes supported by the high-refractive index silicon. Furthermore, 
the attenuation and mode index are significantly altered by conductivity changes in 
the silicon. An amplitude modulator and phase modulator have been proposed using 
these results. Predicted high attenuations in the device may be reduced signifi- 
cantly with a silicon dioxide buffer layer. 


INTRODUCTION 


A need has arisen for direct optical modulation technology. This need has 
arisen from the search for faster digital switches, higher capacity data channels, 
and light, compact data preprocessing equipment for satellites. One promising 
technology that has been examined is the modulation of a guided light wave via 
photoconductivity changes in a semiconductor cladding. 

Computer modeling studies on four- and five-layer, silicon-clad, planar di- 
electric waveguides indicate that the propagation characteristics can be altered 
by changes in the complex permittivity of the silicon and in the thickness of the 
silicon. Using these predictions, an intensity modulator and a phase modulator 
based on photon- induced conductivity changes in the semiconductor cladding have 
been studied. 


SEMICONDUCTOR-CLAD WAVEGUIDES 


The four-layer planar waveguide structure under consideration is shown in 
Figure 1, where it is assumed that light is propagating in the dielectric (N^) 
and all materials are lossless except for the semiconductor (N^) • We desire to 
solve for the complex mode propagation constant (a + j 3) . 


^’^Research sponsored by NASA-Langley Research Center under Grant NSG-1567. 


217 



One technique ( 1 ) extends Maxwell’s equations and boundary conditions to 
numerically solve a transcendental equation relating the attenuation constant 
(a) and phase constant (B) to the material types and thicknesses of the wave- 
guide structure (hereafter refered to as PROGRAM WAVES). 

A more efficient method ( 2 ) utilizes a matrix representation of Maxwell’s 
equations, field solutions and boundary conditions in each waveguide layer (here- 
after refered to as PROGRAM MODEIG) . The matrices are multiplied and a character- 
istic matrix for the entire structure is obtained which yields the attenuation 
constant and phase constant. 

The waveguide consists of a semi-infinite glass substrate, a polystyrene core 
of thickness 1 micrometer, a silicon cladding of .01 micrometer to 10 micrometers 
in thickness, and a semi- inf inite layer of air. Each material is characterized by 
a complex relative permittivity, e; a free space wavelength of 632.8 nanometers is 
assumed and material parameters are shown for this wavelength (Fig. 1 ) Layers 
N^, and N4 are lossless dielectrics, so is real; however, at optical frequencies 
the permittivity of the silicon (N2) is complex ~ ^*r complex 

part is a linear function of the conductivity (^’’^ ~ o/b^e^) . 

PREDICTED CHARACTERISTICS 


The curves presented in Figures 2 and 3 were generated by repeated use of our 
own PROGRAM WAVES and later confirmed with PROGRAM MODEIG. The silicon cladding 
was varied from .01 micrometer to 10 micrometers and the complex mode propagation 
constant was calculated. The expected result was that as the cladding thickness 
was reduced to zero, the attenuation decreases to zero in a well-behaved manner; 
however, the results were not well-behaved when the silicon thickness falls below 
1 micrometer. The curves are similar to exponentially damped sinusoids. Extreme 
3/K variations correspond to median values in the a-curve, and extreme a variations 
correspond to median B/K values. By increasing the conductivity of the silicon 
cladding, the amplitude of the curve oscillations decreases slightly, and the a- 
curve shifts vertically to a higher attenuation. The percent change in attenuation 
compared to dark conditions (0 = for different conductivities (a = 1.1 a^, a = 

1.25 0Q, a = 1.5 Oq) is shown in Figure 4 . The B/K curves shift for a conductivity 
change, as well, although not in such a well-defined manner. The percent phase 
shift compared to dark conditions (o = a^) for different conductivity changes is 
shown in Figure 5 . These effects will be used for intensity and phase modulation 
in a device where a signal source induces photoconductivity changes in a thin sili- 
con cladding on a waveguide, thus modulating a coherent beam in the guide. 


FIELD CALCULATIONS 


We now consider the problem of correlating the local maximum /minimum points 
on the attenuation and mode index curves with the electric and magnetic field 
distributions in the waveguide layers. Results Indicate that the presence of a 
thin silicon film (< 98 a; has little effect on the wave function profiles; the 
profiles are similar to those of the three-layer lossless structure (air-dielectric- 
substrate) . For the thick silicon film structure, however, the lowest order mode 
of the lossless three-layer structure couples to the modes associated with the 
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semiconductor film (the high-refractive index silicon behaves as a waveguide) . 
Furthermore, the coupling between the modes supported by the three-layer lossless 
structure and the high loss TE’ modes of the silicon waveguide determines the 
attenuation and phase of the complete four-layer structure. 

Our results can be described as periodic coupling between the guided mode and 
other leaky modes of the same guide. First, we examine the partial structure con- 
sisting of a silicon guiding region surrounded by semi-infinite layers of air and 
polystyrene. The attenuation and mode index are shown in Figures 6 and 7. We note 
a phase match condition between the modes of the partial structure (air, silicon 
guide, polystyrene) and the TEq mode of the complete waveguide at cutoff thicknesses** 
for successively higher order modes of the partial structure. The sharp peaks on 
the attenuation curve, for the four-layer structure, occur whenever the guided wave 
is strongly coupled into the high-loss modes of the silicon partial structure; 
conversely, the sharp nulls of apparently zero coupling efficiency occur at thick- 
nesses midway between the values for two adjacent leaky modes of the partial struc- 
ture. This is similar to the results of power transfer calculations for linearly 
tapered directional couplers (3, 4). The abrupt transitions on the mode index 
curve of the complete structure occurs when the phase match condition is satisfied 
and the guided wave couples into successively higher order modes of the partial 
structure. 

We now consider our results in terms of the electric and magnetic field distri- 
butions at the local maximum /minimum points on the attenuation vs. silicon thickness 
curve. The real part of the TEq mode electric field profile in the transverse direc- 
tion is shown for a cladding thickness, t 2 = .007 micrometers (Figure 8(a)). We recall 
that this thickness is below the cutoff value for the silicon waveguide structure; 
note that the x^ave function profile is not appreciably distorted. 

For the first local minimum (t 2 = .05ym), we observe that the field strength at 
the silicon-dielectric interface approaches zero. We also note the exponentially 
decaying solutions in the outer, semi-infinite layers as expected. For the first 
local maximum (t 2 = .09ym), we observe a sharp peak in the wave function profile at 
the silicon-dielectric interface as coupling to the TE’j^ mode of the silicon guide 
occurs (Figures 8(b) and 8(c)). 

Similar behavior is noted for the next local minimum/maximum pair, Figures 9(a) 
and 9(b), (t^ = .13ym and t^ = -18yra). The field strength is effectively zero at the 
silicon-dielectric interface for the local minimum, and a sharp peak in the wave 
function is evident for the local maximum as coupling to the TE ’2 iriode of the silicon 
guide occurs. We also note that the field begins to oscillate in the semiconductor 
cladding as we couple into the higher order modes of the silicon waveguide structure. 

Again, for the next local minimum/maximum pair (t 2 = .22ym and t 2 = .26ym) the 
number of field oscillations in the silicon increases. (See figures 9(c) and 9(d).) 
The increase or reduction in field strength at the silicon-dielectric interface is 
also apparent. At this local maximum (t2 = .26ym), there is the sharp field peak 
at the interface, but the field decays rapidly through the dielectric indicating 
almost complete energy transfer. For the other maxima cases considered, the sharp 

* TE’. denotes guided modes in the semiconductor and TE^ denotes guided modes in 
the dielectric. 

** Cutoff for the silicon guide occurs when (3/K)c,« ^ ^ 

^ polystyrene 
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peak was evident at the interface, however, a sizable field was still present in the 
polystyrene dielectric. This indicates that the local Tnaximum (and likely, minimum) 
values used for the calculations are not the precise values as in the former case 
(Figures 9(a) and 9(b)). 

The field plots indicate, then, that the attenuation and mode index of the 
four-layer structure may be explained quite simply as a coupling between the basic 
three-layer lossless waveguide (air-dielectric-substrate) modes and the high loss 
TE’ modes of the silicon guide. For a local minimum on the attenuation-thickness 
curve, the field at the semiconductor interface is zero, and for a local maximum, 
a sizable field is set up at the interface which decays rapidly through the guide 
and substrate. Finally, the number of field oscillations in the silicon increases 
as we couple into the higher order modes of the partial structure. 


REDUCTION IN ATTENUATION THROUGH USE OF FIVE-LAYER 

STRUCTURE 


Thin dielectric buffer layers have been used to lower the attenuation losses 
of metal-clad dielectric waveguides (5). These layers are placed between the di- 
electric core and the metal, and act as buffers to remove a large portion of the 
field from the metal claddings. We now consider the effect of an Si02 buffer layer 
on the attenuation vs. silicon thickness characteristics. 

The result for an Si02 buffer layer (e^ = 2.12) of several different thick- 
nesses is shown in Figure 10. We note the familiar damped sinusoidal behavior and 
the corresponding reduction in attenuation. 

o 

The result for an Si02 buffer layer (tgj_Q2 = 2000A) of several different per- 
mittivities is shown in Figure 11. Again we note the damped sinusoidal behavior 
and the corresponding reduction in attenuation. 

Our studies indicate, then, that the attenuation may be reduced significantly 
with an Si 02 buffer layer while still preserving the oscillatory behavior of the 
attenuation curve; more effective reduction is accomplished with a lower permittiv- 
ity buffer layer. Also, a buffer layer increases the 3/K values slightl}?' but 
decreases the amplitude of the oscillations on the 3/K- thickness curve. 


CONCLUSIONS 


Computer modeling studies on four-layer silicon-clad dielectric waveguides 
indicate that the attenuation (a) and mode index (3/K) behave as exponentially 
damped sinusoids as the silicon thickness is increased. The observed effect can be 
explained quite simply as a periodic coupling between the guided modes of the loss- 
less structure and the lossy modes supported by the high-refractive index silicon. 
Furthermore, the attenuation and mode index are significantly altered by conducti- 
vity changes in the silicon; an amplitude modulator and an intensity modulator have 
been proposed using these results. Predicted high attenuations in the device may 
be reduced significantly with a silicon dioxide buffer layer between the semiconduc- 
tor and the polystyrene guide. 
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Experimental confirmation of the predicted characteristics is still necessary. 
A number of thin-silicon film waveguides have been RF sputtered but attenuation 
measurements to verify the damped oscillatory behavior are forthcoming. Conduc- 
tivity variations of the silicon should demonstrate the modulation capabilities. 
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Figure S.- Wave function profile 
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Figure 10.- SiO„ buffer layer (e = 2.12) of different thicknesses 
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